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Abstract 

In the present paper we shall study (2 + 1) dimensional Zjsi gauge theories on a lattice. It is shown 
that the gauge theories have two phases, one is a Higgs phase and the other is a confinement phase. 
We investigate low-energy excitation modes in the Higgs phase and clarify relationship between the 
Zn gauge theories and Kitacv's model for quantum memory and quantum computations. Then 
we study effects of random gauge couplings(RGC) which are identified with noise and errors in 
quantum computations by Kitaev's model. By using a duality transformation, it is shown that time- 
independent RGC give no significant effects on the phase structure and the stability of quantum 
memory and computations. Then by using the replica methods, we study Zn gauge theories with 
time-dependent RGC and show that nontrivial phase transitions occur by the RGC. 

1 Introduction 

In the last few years, discrete gauge theories have got renewed interests as a possible device for 
the quantum computations, a quantum computer. This idea was first proposed by Kitaev in his 
seminal paper pjj, and after that there appeared interesting works on this idea|21 El E E] . One of 
the most difficult problem of making a quantum computer and performing quantum computations 
fault-tolerantly is the stability of the quantum states which participate in quantum memory and 
computations. There must be a (large) energy gap between these states and others in the system 
and also mixings of these states must be suppressed by certain effects or selection rules. Then 
one can conceive that topological interactions such as the Aharonov-Bohm(AB) effect may play 
an important role there. The AB effect in the two spatial dimensions gives nontrivial statistics to 
particles with gauge interactions, i.e., anyons. The groundstates of the anyons are degenerate if 
the space is a torus and almost no mixing occurs between them because of the topological quantum 
number. Whereas the gauge symmetry should be descrete in order to avoid long-range interactions 
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besides the topological ones. Kitaev's model is based on the anyonic excitations in qubits system. 
However its detailed relationship to the gauge theory is still missing. 

In this paper, we shall study discrete Zjv gauge theories in (2 + 1) dimensions. There are two 
phases in these model, one is a confinement phase and the other is a Higgs phase. We show that 
Kitaev's system of qubits corresponds to some limit of the Higgs phase of the Z2 gauge models. 
Stability of Kitaev's model is closely related with the confinement-Higgs phase transition of the Zjq 
gauge models. 

This paper is organized as follows. In Sec. 2, we study Z^q gauge theories which appear as a result 
of spontaneous breakdown of C/(l) gauge symmetry. We clarify the relationship between the gauge 
system and Kitaev's model for quantum memory and computations. In Sec. 3, low-energy excitations 
in the Higgs phase are investigated. There appear anyonic excitations, magnetic vortices and dyons 
in a natural way as in the spontaneously broken gauge systems in the continuum space 6 . In Sec. 4, 
phase structure and effects of the (static) random gauge couplings(RGC) are investigated by using 
a duality transformation. The Z^q gauge systems are transformed to spin systems which are more 
tractable than the gauge systems. In Sec. 5, effects of the time-dependent RGC are studied by the 
replica methods. It is found that nontrivial phase transitions occur as the RGC varies. Section 6 is 
devoted to conclusion. 



2 UiX) gauge theories 

Let us start with the following UiX) Abelian gauge-Higgs model on a 2-dimensional(2Z?) square 
lattice. Hamiltonian is given by, 

link plaquette x link link 

+M J2 4^-. - 7 E ^l+.t^«V- + E + (1) 

X link X 

where Uxi is the L'^(l) gauge field on the link {x,i){x = site, i = 1 or 2) and E^i is the conjugate 
electric field. The Higgs field (j)x G C/(l) carries U{1) charge N whereas the charge of the fermion 
tield^ipxifx) is q{—q) which is an integer. H!^ is the conjugate field of (j)x, the gauge coupling is g and 
the fermion mass is M. Other notations are standard. We are interested in the case N ^ 1. In this 
case there are two phases in the model, one is the Higgs phase and the other is the confinement phase. 
In particular in the limit — > 0, the gauge field Uxi is restricted to the pure-gauge configuration 
and the model reduces to a Hamiltonian description of the classical 31? XY spin model plus the free 

^We often call ipx and ipx fermion because they satisfy fermionic anticommutation relations. As a result of the 
gauge interactions, they obey anyonic statistics in the Higgs phase. See later discussion. 
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fcrmion system. The classical 3Z) XY model exhibits a phase transition from the magnetized phase 
to the disordered phase at a critical coupling Kc- On the other hand for large k, quantum fluctuations 
of (j)x are suppressed and low-energy excitations of the gauge and Higgs fields are restricted as 
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1, 



up to (time-independent) local gauge transfomation. Then we can put 

where the gauge operator Z^i is explicitly given as follows by {N x N) matrix, 

^ I Q ■ ■ \ 

eifi • 



(2) 



(3) 







2x(]V-l)., 

e « 



(4) 



Corresponding to the above representation of Z^i^ we introduce "conjugate matrix" X^i as follows, 

'^010 • \ 
10 



0-01 

1 • • 



(5) 



One can easily verify the following commutation relations, 

XxiZxi ^ e~'^ZxiXxi, XxiZyj = ZyjXxi ioy {x,i) ^ {y,j). (6) 
The electric term in Eq.(^ is reduced to the following term in the reduced Zpf gauge-field space. 

The above result can be shown by using the eigenstates of the electric fields as basis vectors. Let us 
define an "empty state" |0) as 

m = 0, (7) 

where we have omitted link index for notational simplicity. By using the following commutation 
relation, 

[E, U] = f/, (8) 
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we can show 

E^U\0) ^U\0) = \1). (9) 

Then the gauge field U is the rising operator of the electric field. 

For the Zn case, we also define "empty" state for the X operator, 

^|0)x = |0)x, X^O)x = \0)x. (10) 

The state \Q)x can be expressed by the eigenstates of the Z operator \k)z, Z\k)z ~ e'^^''\k)z, 

\0)x^^(\l)z + \2)z + --- + \N)z). (11) 



Then one can easily show, 

Z\Q)x = \l)z + e^'\2)z + • • • + e2-^'|7V)z, (12) 

and therefore 

XZ|0)x = e^*Z|0)x, (13) 
(X + Xt)z|0>x = 2cos(^) Z|0)x. (14) 

From Eas. Hl()(l and H14(l . Z is the lowering operator of X + X'^ and therefore, 

i?2,^_(X„ + X^), (15) 

up to irrelevant additive and multiplicative constants. Then for large k, the C/(l) gauge theory 
reduces to the following Zx gauge theory, 

Rt = Hz + Ht + H'^z^ 

Hz = -Ai^X„~A2^^ZZZ + H.c., 

= -7E^i+«^"V + ^^E^i^-+H-C., (16) 

where Ai and A2 are coupling constants of the Zx gauge theory and they relate to the U{1) gauge 
couphng g'^ as Ai ~ .g^ and A2 ^ 1/5^. 

The above "derivation" of the Zx gauge theory (|16|l from the U{\) gauge system is rather 
sketchy but it might be useful for realization of discrete gauge systems in architecture of the quantum 
computers. For example, spontaneous breaking of U{1) gauge symmetry occurs in the supercon- 
ductivity. In most of the superconductors including the high-temperature ones, the "Cooper pair" 
carries electric charge 2e. Then a discrete Z2 gauge system close to the present one might be realized 
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in some superconductors. The Hamiltonian H16(l is directly obtained from the path- integral formal- 
ism of the Zn gauge theory on 3-D lattice by taking the continuum limit of the time-like direction. 
In the 3-D Zn gauge theory, there exist two phases, i.e., confinement and Higgs phases as we show 
later on. Phase transition occurs at a certain critical coupling (Ai/A2)c- In the original U{1) gauge 
theory, there exists a critical line connecting the XY phase transition at (g — 0, k — Hc) and Zn 
gauge phase transition at {g ~ gc, n ~ oo) (see Fig.l)[7]. 

Physical state of the system (I16II must be gauge-invariant and this condition is given as follows, 

= ( n Xy,)e-"-^'^^^''^--^'^'^-\ G,\phys) = \phys), (17) 

{y,i)ex 

where {y, i) x denotes 4 links emanating from site x and Xyi — Xyi for y = x whereas Xyi — 
for y — i — X. It is proved that Gx is the gauge-transformation operator at site x and the Hamiltonian 
+ H'^ in Eq. commutes with Gx ■ 

Recently Kitaev proposed a 2-dimensional qubits system for fault-tolerant quantum memory 
and computations pp. This system is closely related to the Z2 gauge theory and contains "anyonic 
excitations" . The system is defined on a torus and the Hamiltonian is given as follows in our notation, 

Hk = -Y, n Xy,-J2zZZZ, (18) 

x {y,i)ex pi 

where Zxi and Xxi are explicitly given by the Pauli matrices in the Z2 case, i.e., Z = and X — a'^ . 
The groundstates and excited states of the Hamiltonian (|18|) are easily obtained since the first and 
second terms of (|18|l commute with each other. The groundstates satisfy 

Yl Xy,\GS)K = \GS)k, Y[z\GS)k = \GS)k, (19) 

{y,i)Gx pi 

for all sites and plaquettes. The groundstates are four-fold degenerate on the torus, as we explain 
in the following section. These degenerate lowest-energy states form basis for quantum memory T. 
The first excited states are explicitly given by 

Jl Xy,\lst)K^-\lst)K, or ]JZ\lst)K^-\lst)K, (20) 

{y,i)Gx pi 

for some specific site or plaquette and otherwise they satisfy Ea. (|19|l . It is not so difficult to see that 
Kitaev's model is equivalent to the model H16|l with N = 2,j — 0,M — 2,q— 1 and Ai = 0, A2 = 1. 
With these parameters and the physical state condition H17|l . the groundstates of the gauge model 
are given as 

l[z\GS)z = \GS)z, i^ii^x\GS)z^O, 'fl^x\GS)z ^ 0, (21) 
pi 
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for all plaquettes and sites. From (|17|l . the second and third conditions of H21|l mean Yl(^y i)^^ ^yi\GS)z — 
\GS)z- On the other hand, the first excited states of the gauge system are given by, 

n^llst)^ = -\lst)l, or = \lst)t or = \lst)%, (22) 

pi 

for some specific plaquette or site. From Ea. (|17|l . the second condition in H22|l is equivalent to 
Y\^y i)(zxXyi\lst)^ — — |lst)^ for e'^^^^^^''\lst)^ — — |lsi)^ and energy increases by 2 because of the 
mass term in (116(1 with M = 2. Similarly for the other fermion ip^. In the original paper by KitaevQ, 
relationship between his model and gauge theories was slightly discussed but full relationship was 
missing. In the following sections we shall study phase structure of the present Zn gauge model, low- 
energy excitations, effects of random gauge couplings, etc. All these discussions give an important 
insight to the stability problem of Kitaev's model. 

3 Low-energy excitations in the Higgs phase 

As we show in the following section, there are two phases in the Zjv gauge theory in (2+1) dimensions 
Hz in p6(l . For large A2/A1, fiuctuation of the gauge field Z^i is small and the Higgs phase is realized 
whereas for small A2/A1, the gauge field Z^i fluctuates strongly and the confinement phase is realized. 
The Higgs phase of the model can be used for a quantum memory. Coupling of the matter fields ipx 
etc. enhances the Higgs phase. 

Let us study the model on the torus and focus on the Higgs phase for large A2/A1. In particular 
for Ai — 0, the groundstates are given by Ea. l21|l and low-energy excited states are particle states of 
ipx, 'fix and states of plaquette magnetic excitation or vortex, i.e., Yipi Zllst)"^ ~ — |lst)^ for specific 
plaquette. As we study the model on the torus, we have the following "trivial" identities 

n n =1 n n^^^ (23) 

all sites (a,i)ea:: all pi's pi 

and therefore the above excitations must appear in pairs. As the groundstates satisfy Ea. (|21|l . there 
is no mangetic flux in each plaquette. Then one may think that the groundstate is unique. However 
this is not the case. There are two nontrivial cycles on the torus, and let us call them a-cycle and 
b-cycle, i.e., noncontractible closed paths. We introduce the dual lattice in the usual way, and choose 
certain noncontractible closed loops on the original and dual lattices. We use notations such that 
Cf (C^) for a suitably chosen closed loop corresponding to the a-cycle(b-cycle) on the original lattice 
and C^(C^) for a loop corresponding to the a-cycle(b-cycle) on the dual lattice. Later discussion 
does not depend on the choice of the loops. Then we define the following operators, Za, Zi,,Xa and 
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^a=n^-' ^6=n^-' (24) 

where XxiS in cross and similarly for X^. These operators are obviously invariant under 
gauge transformation and commute with Hz when Ai — 0. Furthermore they satisfy the following 
commutation relations, 

XaZb^e'^'ZbXa, XbZa^e'^'ZaXb, (25) 

and otherwise commute. Therefore the groundstates are eigenstate of the the operators, e.g., Za 
and Zfc and they are iV^-fold degenerate. This result holds even in the presence of the fermions ijjx 
and ipx since Za and Z^ commute with H-j for vanishing Ai in Ea. ljlfill . 

Fermions f/'a: and ipx move in an unfluctuating "background" field of ZxiS with vanishing magnetic 
field. However they distinguish the above A^-^-fold degenerate Z's groundstates. In fact while iI)x{oy 
Lpx) fermion moves along a closed loop of the a-cycle, it acquires phase factor which is an eigenvalue 
of {Za)'', and similarly for the b-cycle. Then the Higgs phase is a "topologically ordered" phase. 
The iV^ groundstates work as qudit for quantum memory and the quantum states of the qudit are 
distinguishable by using matter fields like ipx- 

Let us discuss excitations in detail. As we explained above, the fermions must appear in a pair. 
Two-fermion state at sites x and y is explicitly given as, 

\F;Cxy)^4{l[z^)^l\GS)z, (26) 

where Cxy is a certain path on the original lattice connecting x and y, and the state (|26|) obviously 
satisfies the physical-state condition l|17|l . On the other hand two- vortex state at dual sites x* and 
y* is given as, 

\V;Cx^y^)^[ n x)\GS)z, (27) 

where Cx'y is a certain path on the duanattice connecting x* and y* and X's in H27|) are on the links 
crossing Cx*y* (see Fig. 2). This state is also a physical state. Other physical excitations are produced 
by appling the gauge-invariant operators in Eas. H26|l and H27(l succsessively on the groundstates. 

Fermionic excitations and magnetic vortices satisfy a nontrivial statistics. This is an Aharonov- 
Bohm effect of the Z^ gauge theory. To sec this, we consider the state like 

<{ n n x)\GS)z, (28) 

^2^2 
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and assume that the paths Cx-^yi and Cx'y* do not entangle with each other. Let us move the ip^ 
fermion at xi around the vortex at X2 once counterclockwise (and not j/2) ^-nd then return it to the 
original position xi. The resultant path C'^^y^ encircles x^ once and C'^^y^ and Cx^y* cross with each 
other. Then the state can be written as 

<i n n x)\GS)z=ii( n n n ^)\gs)z, (29) 



^'^ivi '-^^5»2 ^"'ivi '-^closed '-^'=2^2 

where Cdosed is the closed path (C^^j^^ • C~^yJ which encircles x^ once and has a single common 
link (or odd number of links) with Cx*y»- Because of the nontrivial commutation relation between 
Zxi and Xxi and l|21(l . the resultant state differs from the original one by the phase factor e""^*, 

n n n ^)iG5)z=e^v^,( n n x)\gs)z. m 

^'closed C^'yj '^^ivi C^Jb* 

The above anyonic properties of the low-energy excitations are closely related with the ground- 
state degeneracy. In the continuum spacetime, a Chern-Simons(CS) gauge theory is often employed 
for describing anyons which are a nontrivial representation of the braid group. In anyon systems 
on a torus, movement of an anyon along noncontractible loops like the a-cycle and/or b-cycle is 
a nontrivial element of the braid group. On the torus, the zero modes of the CS gauge field play 
an important role and the groundstate wave function of anyons becomes multi-component because 
of the zero modes[H]. Similar phenomenon occurs in the present Z^r-gauge system as we explained 
above. 

One may conceive that the system has dyonic excitations as in the continuum theories 6 . The 
answer is positive. Dyon dx with "electric charge" Qe and "magnetic charge" R is described by the 
following Hamiltonian, 

Hd^-J2 ^I^. + H.c, (31) 

where we assume that the fields dx and dj, themselves satisfy the fermionic commutation relations 
for simplicity. The link (xi) is associated with the link (xi) and defined as follows, 

f (x + 1,2) for i = 1, 
link (xi) = { (32) 
[ {x + 2, 1) for i = 2. 

From the above definition H31() and H32|l . it is obvious that the electric charge of the dyon da;is 
located at the site x whereas its magnetic charge R is located at the nearest-neighbor plaquette (see 
Fig. 3). Regularization is naturally introduced by the spatial lattice. It is not so difhcult to show that 
the above dyon satisfies nontrivial representation of the braid group and there appears the phase 
factor like — exp(± '^iQE+R)-^i -j .^^jj^gj^ f^^^Q dyons interchange with each other. 
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When we turn on the parameter Ai in Hz the operator Za and Zjj do not commute with 

Hq- anymore and therefore degeneracy of the groundstate disappears. This stems from the fact that 
because of the term Ai "^X^i, Z^i becomes dynamical and it fluctuates quantum mechanicaUy and 
then genuine anyonic properties of the low-energy excitations break down. However for small Ai, 
there is still an energy gap between the iV^ "groundstates" with fine structure and the other excited 
states. Furthermore, these N"^ states are far apart with each other in the quantum-mechanical 
configuration space and are hardly mixed if the torus is sufficiently large. Therefore the system 
with small value of Ai is still suited for a quantum memory as Kitaev suggested first. However 
as Ai increases, a phase transition occurs as we show in the following section. In the new phase, 
a confinement phase, the gauge field fluctuates randomly and the system is useless as a quantum 
memory.^ 

4 Duality transformation, phase transition and random gauge 
couplings 

In the previous section, we discussed that for the quantum memory and commputations the Higgs 
phase must be realized in the present system. In this section we shall study the phase structure of 
the gauge-theory model Hz in (|16() . To this end, the system is deflned on a large spatial square 
lattice. We shall perform a duality transformation which transforms the gauge-theory model into 
a more tractable spin model. For the Z2 gauge theory, the duality transfomation is discussed in 
Kogut's review article [TU). 

Let us consider the pure gauge system Hz in H16|l with the physical state condition. 

By solving the above condition 1)^^3(1 . the operator Xx2 is given as follows by the remaining operators, 

Xx2 = xl^xl _^xl_2^iXl_2_i ■■■ . (34) 

As the "conjugate" operators Zx2 of Xx2 commute with the Hamiltonian Hz, we can set it as a 
constant, Zx2 = 1- 

Then we introduce the following dual operators Wx* and T^. which reside on sites of the dual 

^It is very interesting to see that similar gauge-theory argument can be appHed to neural network models for 
brain|2]. There Higgs phase corresponds to good brains and the confinement phase to dementia. 
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lattice, 

= n ^V^^ 

V,, ^Y[x,_2i,i, (35) 

where (yi) G x* denotes 4 links to the plaqutte on the original lattice which is dual to the site x* of 
the dual lattice. From the definition (|35|l . one can easily verify relations like, 

V,.WI, =e-'^wl.V,., (36) 



K-KI*-2 yi.l4*-i =^.2. (37) 

From Eas. H3t)|) and H37|) . the Hamiltonian Hz in pt)|) can be rewritten in terms of Vx* and Wx-, 

Hz = -Xi J2 14.'^^i._,-A2^W^.. +H.C. (38) 

£c*,i=l,2 X* 

The above quantum Hamiltonian (|38|l is nothing but that of the 3D classical Zn Ising model(the 
clock model) which is obtained by the transfer-matrix methods and taking the continuum limit of 
one direction. 

The Hamiltonian (|38|l is more tractable than the original one (|16|l . There are two phases, i.e., 
ordered and disordered phases, and a phase transition occurs as the value A1/A2 varies. For small 
A1/A2 limit, the groundstate is given by 



w^,.|o)5 = w^i- 10)5 = 10)5. 



(39) 



In the representation. 



1 

e^' 



e " 



(40) 
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the above groundstate \0)s is explicitly given as, 







(41) 



vv 

For small but nonvanishing A1/A2, the groundstate is obtained by the usual perturbative calculation, 
and the term Vx'V^._^ tilts nearest-neighbor Wx* and Wx*-i by ±7^, respectively. In this phase, 
there is no "magnetization", i.e., 

s(0|F..|0)s = 0. (42) 

Low-energy excitations are given by, 

1 ^ 1 



Vx^\0).i 



e'^'->i.|0)s, 



(43) 



where 2-vector p is a momentum and Ns is the number of the sites. Excited energy of the above 
states (|43ll can be easily calculated and obtained as follows for small A1/A2, 



E = 2X- 



:( — (f))- 



(44) 



From (|44|1 . the energy gap is a decreasing function of N. 

For large A1/A2, on the other hand, the groundstate of the spin system H38() is given by 

V^..|6)5 = e^"|0)s, yj. 16)5 = e-'" 10)5, 



(45) 



where e*" € and therefore it is A'^-fold degenerate.^ There is a nonvanishing magnetization for 
large A1/A2, 

s(6|^..|6)s7^0. (46) 

From Eqs.(|S3), ^ and the discussion in the previous section, it is obvious that vortex condensa- 
tion occurs in the gauge-system state corresponding to \0)s- This means that for large A1/A2 the 
confinement phase is realized and therefore the gauge system does not work as a quantum memory. 
This result is important for the architecture of the quantum computer. 

It is interesting and also important to study another type of disturbance for realization of the 
Higgs phase, or, a good quantum memory, i.e., the effect of random gauge couplings which corre- 
sponds to noise and errors in quantum computations. In this section, we consider static random 

^This result of the spin system does not mean that the original gauge system has the degenerate groundstates. 
Actually from I16i . the groundstate satisfies X^ilGS; gauge) = 1 for all links (xi) in the large A1/A2 limit. 
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gauge coupling(RGC) with random variables Tx* G Zn, and the Hamiltonian is given by, 

link pi 

We assume a simple local correlation for the random variables, 

[Tx-T-y.] CX (5a;.y., (48) 

where [• • • ] denotes the ensemble average. 

It seems rather difBcult to study the above random gauge system (|47|) . However by using the 
duality transformation H35|l . we can rewrite the Hamiltonian as in the nonrandom case which 
we studied above, 

Hz = ~\i J2 K*V;t_,-A2^T,.W^,. +H.C. (49) 

x',i=l,2 X* 

Then we redefine the dual operators Wx' as follows, 

Wx'=Tx'Wx^. (50) 

One can easily verify that the new operator Wx* and the old one Wx- satisfy exactly the same 
operator equations in and also there are no spatial correlations of Wx* because of (|48|l . Then 
is equivalent to the original Hz and the random system has the same phase structure with 
the nonrandom one. The groundstate, excitated states, etc. are different in the two systems but 
there exists one-to-one correspondence between them. This result can be partly seen in the original 
gauge system (|47|l . For static RGC Tx* and Ai = 0, there is a unique Z- field configuration of 
the lowest energy up to local gauge transformations. Vortex excitations are generated by applying 
the string operator (Jl^^^O this lowest-energy configuration as in the nonrandom case. Then 
there is one-to-one correspondence. The perturbative term XiXxi generates a pair of vortices in a 
nearest-neighbor plaquettes from the lowest-energy state as in the nonrandom system. 

In the following section, we shall study random Z2 gauge system with full RGC by the replica 
methods and show that nontrivial phase structure appears. 

5 Replica mean-field theory 

In this section we shall study the d-dimensional random Z2 gauge theories by the replica mean-field 
theory (RMFT). RMFT has been often applied to the random spin systems and spin-glass problems. 
The random gauge theories have been less studied and as far as we know there is no systematic 
studies on the random gauge theories by the replica methods. Numerical Monte-Carlo simulations 
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are also useful to obtain phase diagram in the p — T plane, where p is the concentration of plaquettes 
of "wrong sign" and T is the teniperature |lllll2| . In the previous sections we used the Hamiltonian 
formalism, but in this section we employ the path-integral Lagrangian formalism since the path- 
integral formalism is more suitable for the MFT. 

Let us first study the nonrandom Z2 gauge theory on a d-dimensional lattice by the MFT[T^. 
The partition function Z is given by, 

Z = Tre-'^^, 5" = -^CTcrcrcr- /i^cr^,, (51) 

pi link 

where the Z2 gauge variables a^iH = 1, • ■ ■ ,d) take ±1, Tr means =±ij P is inverse temperature 
and h is an external "magnetic field". It is not so difficult to drive MFT action Sm- To this end, let 
us decompose a^i as axi — Uq + Suxi where C/q is the MF for a^i and 5axi is the fluctuation from 
it jl4j. In terms of the new variables. 



link link 

= -U^Np-U^-2{d-l)Y,ax^ + 2{d-l)U^NL-hY^ax, + 0{{Sax^f), (52) 

link link 

where Np and are the numbers of plaquettes and links of the lattice, respectively and — 
-jzrj^Np. From (|5^ . Sm is obtained as, 

Sm = SU^Np - {2{d - l)U^ + (^3) 

link 

Then it is straightforward to calculate the partition function from Sm in (|53|1 . 

Zmf - Tre-''^- 



-3l3USNp 



2coshp{2{d-l)U^ + h} . (54) 



The "magnetization" m per link is calculated from H54|l as 

1 OZmf 

Tfl ^ 

ZmfNl dh 

= tanh/3{2(d-l)C/o3 + /i}. (55) 
Similarly the free energy is obtained as, 

2cosh/3{2(d - 1)U^ + h}] + WqNp. (56) 



F^-^^iVplog 



Numerical calculation of the free energy F in (|56|l is given in Fig. 4 as a function Uo for vanishing h 
and at various inverse temperatures (3. From Eas. H55(l and (|56|l . it is verified that the magnetization 
m is equal to the value of C/q at stationary points of F = F{Uo). Result in Fig. 4 shows that there 
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is a first-order phase transition as the temperature varies and at low temperature the magnetiztion 
TO is novanishing. 

It is known that there is a second-order phase transition in 3D Zi gauge theory which is dual 
to the 3-D Ising model as we showed in the previous section^UI- This means that the MFT gives 
correct results only at large spatial dimensions as it is well known for the spin systems, etc. However 
we believe that the MFT is still useful for obtaining rough estimations of the physical quantities, 
phase stucture of systems, etc. 

Let us turn to random gauge theories(RGT). We study the random Zi theory with the following 
action, 

Sr = — JpCraaa, (57) 
pi 

where we assume that the RGC Jp has the probability distribution like, 

F(Jp) = -^exp{-^(Jp- Jo)^}, (58) 

with positive parameters J and Jq. We choose ()58(l in which Jp takes continuous real value instead 
of the discrete distribution Jp = ± Jg, because it is more tractable. 

We apply replica tricks to the above RGT and then the partition function is given as, 
„ n 
[n= / {n^^P^(-^p)}Trexp(/3^Jp^n^"+/^^E^")' (^9) 

p p a=l p link 

where a is the replica index which takes a = 1, • • • , n, and we shall take the limit n — > in the final 
stage of the calculation. Because of the replica tricks, integration over Jp can be done for each p, 

We introduce the MP Uoa for cr"j and the glass MF(GMF) Qaf3 for cr"iCrf,. Then the terms in the 
action l|60|l can be rewritten as follows as in the MFT for the nonrandom case H53|l . 

pi p link 

E n n ^ -^^pqu + ^cQip E (61) 

pi p p link 

where we have put C = 

From Eqs.(|Sni), JEOI) and (jHU, 

[Z"] = exp ( - 3/3V27Vp ~ ^JoPNp ^ U^^ + Nl log Tr e^) , (62) 

a</3 a 

where 

L = AP^J^C Qlf3'y"r'^L +pY('^JoCUi^ + h)a^,. (63) 

a</3 a 
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We assume a replica symmetric (RS) solution lor Uoa — Uq and Qap — Q- In the RS case, logTr 
can be evaluated as follows, 

32 72/ 



logTre^ = logTr^i^^:g^|^dzexp(-1^4^z2 + 4/?VW^E^"^ 
-2/32 j2CnQ=' + f}{AJoCU^ + h)Y^ <) 

= logfl + n /" i:)zlog(2cosh/3i/(0))-2?i/32j2CQ3 + O(?i2)V (64) 

J —00 



Dz^dz^-P^, H{z) ^2Jy/CQ^z + AjQCU^ + h. (65) 
V 27r 



where 



From (|64|1 . the free energy is evaluated as 
_pFj, - hmO^ 

ri^O n 

- Np[- ^^{n - l)Q' - iJoPU^ + ^ log Tr 
= Np(^^P^J^Q^-3Jo(3U* + ^ J Dz log{2coshf3H{z))^2(3^J^Q'^y (66) 

The values of MF's Uq and Q are determined by the stationary condition of Fr, 

dFn dFn 

Numerical calculation is necessary for solving Ea. H67() . and the result is given in Fig. 5. 

Let us explain physical meanings of the "order parameters" C/oa and Qap- In the ordinary gauge 
theories with constant gauge coupling, the confinement and deconfinement phases are distinguished 
by the expectation value of the Wilson loop operator W{C)^ 

W{C)^ n (68) 

(xi)ec 

where C is a large closed loop on the original lattice. If the system is in the confinement phase, 
{W{C)) (X e-A'^«=a(C')^ whereas in the deconfinement phase, {W{C)) cx e-Pe"'^^t'='-(C). 

In the RGT, on the other hand, the ensemble average must be taken in order to obtain physical 
quantities. Then order parameter is given by [(VF(C))]. Nonvanishing of the MF [/qq means the 
perimeter law [(W^(C))] oc e"^'^"™''*''''^*^), which indicates that the system is in the deconfinement 
phase, the Higgs phase in the present case. As shown in Fig. 5, the RMFT predicts that the Higgs 
phase exists in the RGT if the fluctuation of the RGC Jp is not so large|15|. Result in Fig. 5 also 
indicates the existence of a "gauge glass" phase^. In this phase, C/qq = whereas Qap 7^ 0. This 
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means [{W{C))] cx g-^''^'^^^) whereas [{W{C)y] cx e-P'="°^^t'='-(C). This prediction of the gauge glass 
is itself interesting but the spatial dimension must be probably large for its realization.^ 

Reliability of the RMFT can be studied as in the usual spin glass models like the Sherrington- 
Kirkpatrick model. Parisi-type solutions for replica-symmetry breaking are also interesting. These 
problems are under study and results will be reported in a future publication. 



6 Conclusion 

In this paper we explicitly showed the relationship between Zn gauge theories and Kitaev's model 
for quantum memory and computations. Zn gauge theories appear as a result of the spontaneous 
breakdown of the U{1) gauge theory with the "Higgs field" of charge N. The Higgs-phase limit of 
the Zn gauge systems corresponds to Kitaev's model. Stability of Kitaev's model was discussed and 
it was shown that the errors or noise represented by the term like ^ X^i induce the phase transition 
to the confinement phase in the gauge-theory terminology. In that phase, quantum memory and 
quantum computations are impossible. Then we studied effects of the RGC which are also regarded 
as noise and errors in quantum computations. Static RGC gives no significant effect on the phase 
structure whereas time-dependent RGC induces phase transitions including that to the gauge-glass 
phase. Application of the present studies to non-Abelian discrete gauge theory is interesting and 
important for quantum computations. 
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Figure 1: Phase diagram of the gauge- Higgs model 
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Figure 2: Path Cx*y* connecting dual sites x* and y* 



18 



X 




x+2,1 



Zx2 'f 




vortex 



charge 



Figure 3: The Hnks (xi) on which operaters. Dyon is composed of electric charge on site 
magnetic vortex on plaquette 




Figure 4: Free energy of the Z2 gauge system by the MFT 
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Figure 5: Piiase diagram of the RGT obtained by the RMFT. The Higgs phase exists at low temper- 
ature(T) and small fluctuation of the RGC. Gauge-glass phase appears at low T and large fluctuation 
of the RGC. In the Higgs phase, the MP's Uq j^O and also Q j^O, whereas in the gauge-glass phase 
[/o = and Q ^ 0. 



20 



